function [N Mean Std Skewness Kurtosis Rho] = Statistic(price_path,num_day)
    temp = log(price_path(2:end)./price_path(1:end-1));
    if num_day == 1
		r = temp;
    else
        r= NaN(floor(length(temp)/num_day),1);
        for i = 1:floor(length(temp)/num_day)
            r(i) = temp(i*num_day);
        end
    end
    N = length(r);
    Mean = mean(r);
    Std = std(r);
    Skewness = skewness(r);
    Kurtosis = kurtosis(r);
    Rho = autocorr(r,5);